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We present a theoretical study of the curvature- induced radiation of surface plasmon polaritons 
propagating around bends at metal-dielectric interfaces. We explain qualitatively how the curvature 
leads to distortion of the phase front, causing the fields to radiate energy away from the metal- 
dielectric interface. We then quantify, both analytically and numerically, radiation losses and energy 
transmission efficiencies of SPPs propagating around bends with varying radii- as well as sign-of- 
curvature. 
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I. INTRODUCTION 

Surface plasmon polaritons (SPPs), coupled modes of 
plasmons and photons, are low-dimensional excitations 
propagating at metal-dielectric interfaces. As such, SPPs 
are confined to lateral dimensions of the order A/ 10, with 
A the vacuum wavelength of light, enabling plasmonic 
devices which are more compact than existing photonic 
equivalents. Owing to their large electromagnetic (EM) 
field intensity at the interfaces, SPPs are highly sensi- 
tive to the surface morphology, thus allowing the real- 
ization of metallic structures capable of controlling and 
manipulating light on the nano-scale. This, in addition 
to advances in nanofabrication technologies, has led to 
a growing interest in realizing ultra-compact plasmon- 
based integrated circuits. Recent studies of plasmonic 
manipulation include SPP waveguiding and bending in 
patterned metallic films [TJ [2], guiding via resonant en- 
ergy transfer in ordered arrays of metal nanoparticles[3], 
as well as demonstrations of SPP prisms, lenses [4 , mir- 
rors, beamsplitters, and interferometers [5]. 

As the need for integration of compact lightwave de- 
vices is growing, it is necessary that we obtain an under- 
standing of the fundamental properties of SPP propaga- 
tion and manipulation in these environments. Of par- 
ticular interest is the development a quantitative the- 
ory of curvature-induced radiative energy loss in SPPs 
propagating at curved metal-dielectric interfaces. This 
enables the determination of SPP propagation efficien- 
cies when the radius of curvature is smaller than or 
comparable to the signal wavelength[6[ [8]. Careful 
analysis of the relations between the propagation effi- 
ciency and the interface curvature is essential when de- 
signing plasmonic devices, as it should set a limit on 
the radius of curvature and, subsequently, on feature 
size in plasmonic-circuits. Several recent studies have 
addressed surface plasmon waveguiding around bends. 
However, these studies are primarily focused on in-plane 
guiding [9j [TQl [TTJ [12] or on using multiple-interface ge- 



ometries such as metal-dielectric- metal waveguides [13] or 
long-range SPP waveguides [14 . To date, SPP guiding on 
a single bent interface remains largely unknown. 

This article presents a theoretical analysis of curvature- 
induced radiation emitted by SPPs propagating around 
bends at metal-dielectric interfaces. In Sec. [Ill we intro- 
duce the geometry of our system, followed by a simple 



qualitative analysis in Sec. [HI] describing the radiative 
nature of SPP propagation around bends. This enables 
us to develop a more rigorous analytical approach for 
calculating losses and propagation efficiencies, which we 
present in Sec. (TV) Section |V| provides a subsequent nu- 
merical analysis, and its results are compared with the 
analytical results. A novel resonator-based method for 



enhancing the transmittance is presented in Sec. VI Con- 



clusions are presented in Sec. VII 
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FIG. 1: (a) The geometry of the curved metal-dielectric in- 
terface. A rounded edge, characterized by a permittivity e*, 
has a bend angle and a finite bend radius R. The bend is 
confined to the region of space shown, with center of curva- 
ture at the origin. The rest of space is occupied by a medium 
with e . Axes x\ and X2 extend along the boundaries between 
regions I and II and regions II and III, respectively. In the 
x\ — X2 plane, regions I and III are semi- infinite. The sys- 
tem is also infinite in extent along the entire z axis. Arrows 
indicate incident and reflected fields in region I, and transmit- 
ted field in region III. In addition, (a) shows the intensity of 
SPPs in greyscale. The intensity distribution is calculated for 
SPPs traveling around a metallic corner (b), using the single- 
mode approximation developed in the text. We also consider 
the complementary configuration, in which SPPs propagate 
around a dielectric void (c). 
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II. GEOMETRY OF THE CURVED 
METAL-DIELECTRIC INTERFACE 

The geometry of our study is shown in Fig. [IJa). The 
system consists of a curved metal-dielectric interface oc- 
cupying region II in space, matching smoothly to planar 
and semi-infinite interfaces in regions I and III. The axes 
x\ and X2 define the boundaries between the three re- 
gions, and the system is infinite in the z direction. The 
rounded edge is characterized by a permittivity e^, with a 
fixed radius R and a finite bend angle 6. The surrounding 
space is characterized by a permittivity e G . Propagating 
SPPs of frequency uj are incident from region I onto the 
boundary at x±, and their counterclockwise transmission 
through region II into region III is analyzed. We initially 
consider the case of SPPs propagation around a metallic 
corner with Re[e^] < and e Q > 0, depicted schemat- 
ically in Fig. [ljb) . This is essentially a segment of an 
infinitely long metal cylinder. We then analyze the com- 
plementary reverse geometry shown in Fig.[TJc), in which 
SPPs propagate around an open dielectric void. 



III. A QUALITATIVE ANALYSIS OF 
BEND-INDUCED RADIATION 



The generalized dispersion relation of SPPs propagat- 
ing at a metal-dielectric (i.e. anisotropic) interface is 



(i) 



where fcy and k± Q ^ are the components of the fc-vector 
parallel and perpendicular to the interface, respectively. 
For a surface-guided mode, we require fen > y/e^co/c 
such that k± Qi i becomes imaginary, hence non-radiative. 
Rewriting this simple expression yields oo/k\\ < c^Je~ Q . 
Simply stated, the phase velocity in the direction parallel 
to the interface, vn = cj/fcy, cannot exceed the speed of 
light, c/y/eZ, in order to sustain non-radiative guiding. 

When propagating around a bend, v\\ acquires a radial 
dependence, with EM fields more distant from the inter- 
face travelling at greater phase velocities. Thus, there 
exists a threshold radius, r*, where the parallel phase 
velocity reaches the speed of light. Beyond r* the EM 
field of the SPP becomes radiative. As a result, SPPs 
can be guided along curved interfaces with negligible ra- 
diation loss as long as the fields are confined near the 
metal-dielectric interface and do not extend beyond r*. 
This claim is verified analytically in the next section. 



IV. TRANSMISSION AND REFLECTION AT 
BENT INTERFACES 

In order to quantify the degree of radiation loss and the 
propagation efficiency of SPPs around the bend we now 
exploit known solutions of Maxwell's equations describ- 
ing angular propagation of EM waves at the surface of 



an infinitely long metal cylinder. These solutions are ap- 
plicable in region II. The magnetic field in the dielectric 
is hence given by 
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where k Q = u^fe^/c and Hn^ is the Hankel function of 
the first kind. The set of mode indices {n}, denoting 
radial excitations, is determined by the metal boundary 
matching equation 
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where ki — LUyfci/c, J n is the Bessel function, and the 
prime denotes differentiation with respect to the argu- 
ment. Since region II comprises only a segment of a full 
cylinder (i.e. < 2tt : ) periodic boundary conditions need 
not be satisfied, and n is therefore not constrained to in- 
teger values. In fact, since we choose the frequency uo 
to be real-valued, one finds n to have a non-vanishing 
imaginary part as well. The latter is a consequence of 
radiation loss and absorption in the bend. 

The azimuthal dependence of the SPP wave is given 
by the standard expression exp[±m0], where (j) is the 
measured from the x\ axis. Thus only solutions with 
Im[n] > are admissible, describing attenuated prop- 
agation. Away from the interface in the dielectric re- 
gion, where k Q r ^> Re[n] is satisfied, we may use the 
asymptotic form of the Hankel function: Hn\k Q r) « 
y/2/iik r exp[i(/c r — (2n + l)7r/4)]. Each mode can now 
be written as 



B^ 
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thus recovering the expected free propagating cylindrical 
wave form. 

In general, Eq. (|3| cannot be solved algebraically, and 
it is necessary to employ a numerical method to find the 
set {n}. One convenient numerical approach employs 
graphical plotting of the right-hand side of Eq. ([3| as 
functions of both Re[n] and Im[n], and identifying its 
zeros in the complex n-plane. 

In principle, the set of solutions denoted by {n} is in- 
finite. In practice, we find that a single mode of this 
set dominates the propagation problem we are analyz- 
ing in this work. The set {n} contains a fundamental 
mode (i.e. the surface plasmon mode), which we label 
with mode index m G {n}. The EM field of this mode 
is concentrated near the metal-dielectric interface. The 
solid line in Fig. [2] shows the dispersion relation of this 
fundamental mode. In the limit of large momentum, it 
asymptotically approaches uo sp 



e , the con- 



ventional limit for surface plasmon modes. Also shown 
in Fig. [2] are the dispersion relations of modes in the set 
{n} with higher radial excitation values (dashed lines). 



3 




10 20 30 40 50 60 70 

Angular momentum number ( Re[n] ) 

FIG. 2: Dispersion relations of various transverse-magnetic 
modes on a metal cylinder with R — 2 /mi. The dispersion 
relation for the fundamental mode (solid line) closely follows 
that of SPPs on a flat metal-dielectric interface (dash-dotted 
line). The horizontal scale of the latter is normalized by the 
factor R, so that the curve is a plot of uj as function of kR, 
instead of k. Dispersion relations of modes with higher radial 
excitations are also shown (dashed lines). These modes are 
less confined to the interface and do not approach co p /^/2 
asymptotically. Higher radial excitation modes with even 
weaker surface confinement exist (not shown) , all lying to the 
right of the three dashed lines shown. 



We see that these modes do not approach uj sp asymptot- 
ically. Compared to the fundamental, these modes are 
also found to be less confined to the interface. As the ra- 
dial excitation number increases even further, the modes 
corresponding to these values of n all lie to the right of 
the three non- fundamental modes shown in Fig. [2] As ex- 
pected from the fundamental mode, Im[m] is the smallest 
of the set {n} since the confinement of this mode to the 
interface is maximal, hence its radiation losses are the 
lowest. This is generally true for metals such as gold and 
silver, where the absorption losses are relatively low [15]. 

As shown in Fig. |2j away from u sp the fundamental 
mode has the smallest angular momentum (« kR) at a 
given frequency. To understand this, we consider the fol- 
lowing: The angular momentum is L = J(r x S/c 2 )<ir 3 , 
where S is the Poynting vector. Compared to other 
modes, the fundamental mode is more strongly confined 
near the surface. Since S for this mode is significant 
only at r c± R, the integrand is minimal (as long as S of 
the fundamental mode is not disproportionately large.) 
However, as the frequency approaches uo spi S increases 
asymptotically, which overcompensates for the smaller 
value of r. This is seen in Fig. [2] as a crossover of the 
fundamental mode and the other depicted modes, such 
that at uo 



its angular momentum is higher than 
that of the less confined modes. 

The radiative nature of the SPP solutions described 
by Eq. ([2| is also consistent with the qualitative analy- 
sis presented in the previous section. This can be veri- 
fied by analyzing the nature of the Hankel function Hn^ . 
The Hankel function is non-oscillatory (depicting non- 
radiative fields) while its argument is smaller than the or- 
der Re[n]. When the argument exceeds Re[n], starts 



to approach its oscillatory (radiating) form, as shown by 
the asymptotic expression leading to Eq. Q. The non- 
oscillatory-to-oscillatory transition occurs when the ar- 
gument and the order are approximately equal. We are 
thus led to conclude that each mode of Eq. Q undergoes 
a non-radiative-to-radiative transition at a radius given 
by k Q r ~ Re[n]. 

The transition radius described above coincides with 
the previously introduced threshold radius, r*. By defi- 
nition, at r = r* we require v\\ = c/ 'y/e^, where the mode 
index n is added to r* since each mode with index n has a 
different threshold radius r* in region II. The phase flow 
in the parallel direction is characterized by exp(zRe[n]0), 
and its associated wavenumber and phase velocity are 
fey = Re[n]/r and v\\ = uo/k\\ = ujr/Ke[n], respectively. 
We thus obtain 



Re[n] 
k Q 



(5) 



verifying that the threshold radius is indeed the non- 
radiative-to-radiating- SPP transition radius. From this 
it follows that the fundamental mode m, with lowest an- 
gular momentum ( for uo < uj sp ,) also has the smallest 
threshold radius for radiation: 



Re[m] 
k Q 



Re 
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The approximation above follows from Re[m] « kR, 
where k = cu^€i€ / (e$ + e Q )/c is the wavenumber of 
SPPs propagating at a fiat metal-dielectric interface. The 
validity of this approximation is confirmed in Fig. [2j 
where it is shown that the dispersion relation of the fun- 
damental mode is very close to that of SPPs at a flat 
interface. 

The above analysis implies that when the majority of 
the SPP field is condfined to a region with r < r^, we 
should expect a negligible radiation loss during propaga- 
tion. On the dielectric side of the interface the SPP field 
decays as exp{— Re[7 G (r — R)]}, with the decay coeffi- 
cient is given by 7 G = uoe ^/— l/Jei + eo)/c. We see that 
the field confinement in the dielectric is therefore charac- 
terized by Re[7 ]. Thus, we conclude that the radiation 
loss is insignificant if G > 1, where G is the exponent 
evaluated at r = r* : 



G (e*, e , ujR/c ) = Re[7 ](r^ - R). 



(7) 



In Fig.[3]we plot Im[m] as function of G. The metal is as- 
sumed lossless (i.e. Im[e^] = 0), therefore Im[m] is related 
solely to radiative loss. We find that Im[m] is a monotoni- 
cally decreasing function of G, as expected from the argu- 
ments above. Moreover, when G = 1 Im[m] = 1.24 ~ 1, 
confirming again that G = 1 is a very good approxi- 
mation of the threshold where SPP propagation changes 
from significantly radiative (Im[m] > 1) to mostly non- 
radiative (Im[m] < 1). 
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FIG. 3: Plot of Im[m] as function of G for uuR/c = 1000 and 
e = 1. The premittivity ei is real and varies from ei — —30 
to -200. 



To quantify the radiation loss and propagation effi- 
ciency around the bend it is now necessary to consider 
the coupling of these metal-cylinder modes to planar- 
interface SPP modes (i.e. propagating modes in regions 
I and III). The magnetic field of the SPP incident from 
Region I is given by 



B = zAexp(z%i — iujt) 



exp[-7 i (xi 
exp[+7 (xi 



R) 
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with A the amplitude and ji = — ujeiy/ — l/(e^ + eo)/c. 
Similar expressions hold for the SPPs reflected into region 
I and the fields transmitted into region III. 

Calculating exact values of the transmission and re- 
flection coefficients is often impractical, since it requires 
matching an infinite number of modes in Eq. ([2| to pla- 
nar SPP modes in regions I and III. The mode matching 
must be carried out along the entire spatial extent of the 
x\ and X2 axes at the boundaries. However, in certain 
cases it is possible to derive relatively simple approxi- 
mate expressions for the required field coefficients. This 
is because although each mode of region II has an EM 
field profile normal to the surface which does not ex- 
actly match the field profile of the incident planar SPP, 
the fundamental mode of region II minimizes this spatial 
mismatch. Non- fundamental modes are characterized by 
EM fields which are less confined to the surface, therefore 
their field profiles deviate more strongly from those of the 
planar SPPs whose fields are always strongly bound at 
the metal-dielectric interface. 

In the short wavelength limit, it is possible to show that 
incident SPP and the fundamental mode have identical 
field profiles near the interface. In this limit, the field dis- 
tribution of the fundamental mode in the radial direction 
can be obtained by taking the appropriate limit of the 
Bessel equation r 2 d 2 f/dr 2 + rdf/dr + {k 2 — m 2 )f = 

to which Jm(kir) and Hm\k Q r) are solutions, respec- 
tively. As ujR/c — > oo, the curvature of the metal sur- 
face becomes insignificant and m approaches kR. With 
r = R+x, in the limit x <C R the Bessel equation reduces 
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FIG. 4: Comparison of intensity profiles of SPPs in region 
I (solid line) and the fundamental mode in region II (dashed 
line), calculated for uoR/c = 800. The profile mismatch is 
barely visible, indicating that SPPs in region I couple pre- 
dominately to the fundamental mode in region II. 



to R 2 d 2 f/dx 2 + (k 2 Q - k 2 )R 2 f = 0, to which the solu- 
tions are exponentials. Thus, noting that / y 2 Q = k 2 — k 2 G , 

Jmikir) ~ exp^r) and Hm\k Q r) ~ exp(— j Q r) which is 
identical to the behavior of the SPP fields in regions I and 
III near the interface. Far from the interface where the 
condition x <C R no longer holds, Jmikir) and H$ (k Q r) 
no longer exhibit exponentially decaying behaviors. How- 
ever, due to the exponential decay near the interface, 
their values are small far from the interface, making the 
mismatch negligible along the entire x\ and x^ axes. This 
is illustrated in Fig. [4j where we show that the intensity 
profiles of the fundamental mode and the planar SPP are 
well matched in the short wavelength limit. 

We are thus led to conclude that in the short wave- 
length limit planar SPPs couple predominately to the 
fundamental mode, and neglect their coupling to all other 
non-fundamental modes. For this reason, under this 
single-mode approximation, it is necessary to consider 
only a small number of modes: the incident and re- 
flected SPPs in region I, the clockwise and counterclock- 
wise propagating fundamental mode in region II, and the 
transmitted SPP in region III. These modes are matched 
at a single point on at each axis, at a distance R from the 
origin, via the standard Maxwell boundary conditions. 
Our analysis above ensures that the boundary conditions 
are then approximately satisfied over the entire extent of 
the axes. The problem of quantifying the propagation 
efficiency has now essentially become one dimensional 
(ID), and it is mathematically analogous to scattering 
from a ID finite potential well 16). However, since the 
allowed m values are always complex, bound-state solu- 
tions in this type of well do not exist. This distinguishes 
SPPs at curved surfaces from waveguide bends enclosed 
on all sides by infinite potential walls [T7]. 

Applying the appropriate boundary conditions to the 
fields at the x\ and x<i boundaries results in expressions 
for the transmittance T and reflectance R: 
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In the presence of significant absorption or radiation loss, 
such that lm[ra]0 ^> 1, these expressions become 
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As uoR/c — > oo, these expressions become exact. 

In general, conservation of energy is expressed by 
1 = T + R + P + A, where P and A are the radiation 
and absorption loss coefficient, respectively. For a metal 
characterized by a real permittivity, the absorption loss 
vanishes, and the radiation loss can be easily calculated 
using the above expressions for T and R. For a lossy 
metal, radiation losses must be calculated independently 
in order to extract the absorption loss from the expres- 
sion above. The radiation loss is obtained by integrating 
the Poynting vector for unit incident flux in region II at 



S • rrdcj). 
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The lower integration limit is set to <po instead of 0, 
since the energy radiated from the surface at <\> = 
propagates at an angle </>o into the far field. Like- 
wise, the upper integration limit is 6 + 0o instead of 
6. In the short- wavelength limit only the amplitude 
of the forward-propagating mode is significant, there- 
fore the radiation losses are well approximated by in- 
tegrating only the counterclockwise propagating mode. 
A stationary phase approximation is used to obtain an 
expression for </>o, using the position-dependent phase 
$ = k Q r + Re[m](j). The change in angle as the wave 
propagates a radial distance Sr is S(j) = Re[m]/(k r 2 )5r, 
giving <j) = Re[m]/(k r 2 )dr = Re[m]/k R. 

Calculations were carried out using typical values of sil- 
ver (ei = — 15 + z0.5) in air (e G = 1) with ujR/c = 800 and 
6 = 90°. Assuming that the metal is lossless (e^ = —15), 
we find that most of the incident SPP energy is trans- 
mitted with T = 0.997, R = 1.19 x 10~ 8 , and P w 0.003. 
When the absorption loss is accounted for the results 
change drastically to T = 0.0516, R = 1.18 x 10" 6 , and 
P « 0.00282. This result implies that absorption is the 
dominant loss mechanism when R ^> A even for met- 
als such as silver with relatively low losses. Surprisingly, 
we find that the overall absorption and radiation losses of 
SPPs propagating at a non-planar interface may be lower 
than the absorption loss of SPPs traveling the equivalent 
arc distance on a flat surface. This counterintuitive re- 
sult comes from the fact that the field inside the metal 
in region II travels an arclength less than OR due to the 
curvature. As a result, SPP fields sample less of the 
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FIG. 5: (a) The upper bound for the transmittance, T n , plot- 
ted for a silver-air interface with bend angle = 90° , as func- 
tion of bend radius R for wavelengths A = 500nm (dashed- 
dotted line), A = 600nm (dashed line), and A = 700nm 
(solid line), (b) Magnified view of the upper bound in the 
diffraction-dominated regime. 



metal volume when propagating on the curved interface 
than when propagating on a flat surface, resulting in the 
reduced absorption. 

In order to evaluate the accuracy of our results, it is 
necessary to quantify the validity of the single-mode ap- 
proximation. As discussed earlier, the single mode ap- 
proximation is only appropriate when the field profile of 
the fundamental mode (represented by the solution to the 
Bessel equation) well approximates the exponentially de- 
caying behavior of the planar SPP. Thus, evaluating the 
mismatch between the Hankel function and the decay- 
ing exponential gives a measure of whether the approxi- 
mation is appropriate or not. We define the normalized 
mismatch as 



Jr 



exp[-7 (r - R)} 



dr 



Jr 



|exp[-7 (r - R)}\ 2 dr 



(14) 



where rj = 0(1). The expression in the numerator quan- 
tifies the field mismatch near the interface. The condition 
A < 1 constitutes a criterion for the validity of our ap- 
proximation. For example, when n = 3, = —15, and 
e Q = 1, A 2 = 0.002 for coR/c = 800, rendering our re- 
sult applicable. On the other hand, for uoR/c = 100 we 
obtain A 2 = 0.3, indicating that the approximation is 
less reliable now, hence the coupling to non-fundamental 
modes can no longer be neglected. 

When the single-mode approximation is not appropri- 
ate, it is still possible to derive a physical quantity from 
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the above analysis. As discussed previously, the mode 
index m associated with the fundamental mode has the 
smallest imaginary part compared to the mode indices of 
non-fundamental modes. Since the wave depends on n as 
exp[±m0], modes with large Im[n] decay rapidly. Thus, 
the transmission in the presence of coupling to nonfun- 
damental modes does not exceed the upper bound of 



exp(— 2lm[m]6). 



(15) 



Here we neglect reflections at the x\ and boundaries, 
thus excluding interference effects. Figure [5|a) is a plot 
of T n as function of R. A peak is clearly visible, mov- 
ing to higher values of R as the wavelength increases. 
To the right of the peak, at large radii of curvature ab- 
sorption losses in the metal dominate, and the maximum 
transmittance decreases with increasing radius. To the 
left of the peak radiation due to the high curvature is 
the dominant loss mechanism, leading to a rapid drop 
in T n . At very high curvature (R < 10/im) there is a 
change in trend, and T n starts to increase with decreas- 
ing R. When calculating the radiation loss per arclength, 
we find that for this range of radii it increases slower than 
elsewhere, allowing T n to increase even as R attains very 
small values. 

The behavior of the transmittance for very high curva- 
ture is a consequence of diffraction of the incident SPPs. 
In general, diffraction is associated with the finite wave- 
length of light, and the diffraction coefficient approaches 
zero in the limit of small wavelength [18]. For SPPs prop- 
agating around bends, a decrease in the radius of cur- 
vature is equivalent to an increase in the effective wave- 
length, X/R. For this reason, SPP diffraction around 
a corner increases as R decreases, resulting in greater 
transmittance for smaller R. Since the diffraction is only 
significant at sharp corners, the transmittance becomes 
weakly dependent on the dispersion of the metal at small 
radii of curvature. Hence, as shown in Fig. |5jb), the up- 
per bound at different wavelengths converges to a single 
value as R approaches zero. 

From the discussion above, it is clear that the non- 
monotonic behavior of the transmittance is a result of 
three competing mechanisms: absorption, radiation, and 
diffraction. Depending on the radius of curvature, one of 
the three mechanisms becomes dominant, creating three 
distinct regimes in Fig. [5|a). However, it is not possi- 
ble to evaluate the potential discrepancies between the 
upper bound model plotted here and the analytical re- 
sults, since our analytical approach is not valid at high 
curvatures. To examine the transmittance for such cases, 
we turn to a numerical study using the finite- difference 
time-domain method (FDTD). 



COMPARISON BETWEEN ANALYTICAL 
RESULTS AND FDTD CALCULATIONS 



We apply here numerical FDTD calculations to study 
the propagation of SPPs about bends, and compare with 
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FIG. 6: (a) FDTD simulation results showing the magnitude 
of the magnetic field for R = A = 630nm. The SPPs are inci- 
dent from the bottom and propagate counterclockwise around 
the bend, (b) Analytical result for the transmittance upper 
bound (black trace) compared to numerical simulations for 
various values of R (red points). The upper bound values 
are clearly higher than the simulation results, confirming the 
consistency of our analytical method. 



our analytical results. Previously, the method of lines 
has been used to study the diffraction and the propa- 
gation of SPPs at a sharp bend with R = 0[19]. Our 
present numerical study examines how the transmission 
efficiency depends on the radius of curvature of the bend. 
The dielectric function in our simulations is given by the 
Debye model 
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With the choice of parameters = 3.90838, e s = 
-25658.4, r = 1.20973 x 10- 14 sec, and a = 1.68931 x 
10 17 sec _1 the Debye model closely matches experimen- 
tally obtained values for silver [20] in the wavelength 
range 400-1 200nm. We implement a nonuniform orthog- 
onal grid, with mesh size in the range A/300-A/40. Dis- 
persive perfectly matched layers [21] are used as absorbing 
boundary conditions throughout our study. The wave- 
length in vacuum is fixed at A = 630nm. A trial simula- 
tion of SPP propagation at a planar metal-air interface 
indicates that the simulated SPPs are well characterized 
by the effective dielectric constant e s = —14.0 + z0.9, 
which is reasonably close to e = — 15.7 + il.l calculated 
using the Debye model. 



A. Ninety-degree bend 

A series of FDTD simulations were performed to an- 
alyze propagation about the 90° rounded edge shown in 
Fig. [6ja). The transmittance, T, is extracted for various 
radii of curvature. As shown in Fig. |6^b), the transmit- 
tance increases with decreasing radius, peaking at R = 
with T « 0.07. Hence, the numerical approach confirms 
the diffraction-dominated small-radius behavior of T pre- 
dicted from our earlier analysis of T n in Fig. [5] Interest- 
ingly, our simulations reveal that the actual transmit- 
tance is reasonably well described by the upper bound 
(calculated now using e s ) even for R < A. Moreover, for 
R > A the discrepancy between the simulation and the 
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analytical upper bound is less than 0.01. This mismatch 
is expected to decrease even further with decreasing cur- 
vature because the coupling to the fundamental mode 
increases as X/R — > 0. We therefore conclude that the 
calculated upper bound is a good estimate for the trans- 
mittance for all R > A. 



B. Bend with negative curvature 

The analytical formalism developed above may also be 
used to analyze the reversed geometry, where the metal 
occupies the outer space, and the SPPs propagate around 
a dielectric void in it as shown in Fig. [ljc). In this 
complementary picture the planar SPP modes are now 
matched to the solutions of a hollow cylindrical void in 
the metal. However, care must be taken when choos- 
ing the appropriate solutions in region II for the mode 
matching. For a dielectric cylinder surrounded by metal, 
the radial solutions in the dielectric are the Bessel func- 
tions. These functions, except for the one confined to 
the interface, have field nodes. Since the existence of 
nodes implies that photons are exchanged between op- 
posing points on the cylindrical interface, these solutions 
are unphysical in an open geometry such as our system. 
Therefore, the only solution represented by the Bessel 
function in region II is the surface confined mode, which 
lacks any nodes. 

As before, the surface confined mode corresponds to 
the fundamental (SPP) mode, which minimizes the field 
mismatch with the incident planar SPP under the single- 
mode approximation. Since the fundamental mode is 
represented by the Bessel function, Eq. ([3| remains un- 
changed and the expressions for the transmittance and 
the reflectance are valid for this reversed geometry. Our 
calculations have shown that this problem is now appro- 
priately analogous to the finite potential barrier model. 
We also find that the absorption loss in region II is now 
greater than the absorption of planar SPPs. From the 
argument above, it also follows that in the single- mode 
approximation SPP propagation around a bend is non- 
radiative. For this reason, SPP transmission through the 
reversed geometry becomes highly efficient when material 
losses and reflections at the boundaries are negligible. 



Our simulations show that efficient propagation is in- 
deed possible. For example, we find that when R = A 
the transmittance is T = 0.73, in sharp contrast to the 
value of T « 0.06 obtained for an equivalent curvature 
in the geometry of Fig. [6^ a). However, when R = 0, 
back-reflection at the bend becomes the dominant loss 
mechanism with R = 0.96, while the radiation loss re- 
mains small as shown in Fig. [7^b). The transmittance 
decreases to T w 0.002. Because of its high reflectance, 
this particular geometry is essentially a SPP mirror and 
may be implemented to construct SPP resonators [22]. 





-2.0 -1.0 0.0 
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FIG. 7: FDTD simulation results showing the magnitude 
of the magnetic field for (a) R = A = 630nm and (b) R = 
when SPPs, incident from the bottom, propagate around 
a circular dielectric void. The radiation loss is significantly 
reduced compared to the propagation around a metallic bend. 



VI. RESONATOR-ENHANCED 
TRANSMISSION AROUND A SHARP BEND 



We have shown that efficient transmission of SPP en- 
ergy around bends is possible when propagating around 
small voids in metals, while SPP propagation about 
metallic bends is highly radiative. In what follows we 
demonstrate a method for reducing radiation losses by in- 
troducing an additional metal-dielectric interface into the 
system, which provides an alternative, low- loss transmis- 
sion channel for SPPs. Consider the geometry shown in 
Fig. |8ja). A cylindrical hole of radius r is placed in close 
proximity to two flat metal interfaces joined by an abrupt 
90° edge. Incident SPPs propagating upward along the 
vertical interface are coupled through their near-fields 
into the void and our of it onto the second interface. 
Positioning of the resonator within the skin depth of the 
metal allows efficient excitation of the resonator modes 
by the incident SPPs, and subsequent outcoupling. The 
role of the dielectric void is similar to that of a dielectric 
microring resonator, and the evanescent wave coupling 
scheme is analogous to optical coupling between micror- 




1 



xOum) 

FIG. 8: (a) Schematic diagram of a SPP microresonator em- 
bedded near a sharp edge. SPPs are incident from the bottom 
along the direction of the red arrow, couple to the resonator 
and couple out again, (b) Numerical FDTD simulation show- 
ing the magnitude of the magnetic field for d = SAnm and 
r = 378nra. A significant portion of the field is seen to couple 
into the resonator. 
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ing resonators and dielectric waveguides [23 . The SPP 
resonator provides an alternative transmission channel 
into which the field couples, leading to a reduction in 
the total SPP energy impinging on the highly radiative 
sharp edge. Hence, the efficiency of propagation increases 
significantly when a resonator is properly incorporated 
into the metal in the vicinity of the bend. Figure [8^b) 
shows the results of a simulation for a typical resonator- 
coupled system. We find that the transmission around 
an infinitely sharp bend in the absence of a resonator, 
(T w 0.07) increases to T = 0.17 when a cavity of ra- 
dius r = 378nm is incorporated at a distance d = SAnm 
from the interfaces. A scattering-theory formalism has 
been previously developed to treat waveguide-resonator 
couplings [24j [25]. However, this method requires ab- 
initio knowledge of the scattering-matrix elements, at- 
tainable by solving Maxwell's equations using FDTD or 
other numerical methods. When applied to the geome- 
tries of our system, we find that this approach is less 
straightforward than our demonstrated method of ex- 
tracting the transmission efficiency directly from FDTD 
simulations of resonator-coupled interfaces. The trans- 
mitted signal levels may be then optimized by fine-tuning 
additional parameters of the system, such as resonator 
shape, the number of (cascaded) resonators and their 
relative positioning near the bend. In addition, incor- 
porating a gain medium[26j [27] (i.e. SPP amplifier) in- 
side the resonator cavity may further enhance transmis- 
sion efficiencies. Future work will address the enhanced 
transmission via SPP resonators in detail. 



VII. CONCLUSIONS 

In summary, we have studied the propagation of SPPs 
at a curved metal-dielectric interface both analytically 



and numerically. We have presented and confirmed a 
physically intuitive picture of the bend-induced radia- 
tion, and how the radiation loss is related to the SPP 
confinement at the interface. In the short wavelength 
limit, we have shown that calculating the propagation ef- 
ficiency is analogous to a quantum mechanical ID finite 
square well which is readily solved to obtain the expres- 
sions for energy transmission and reflection coefficients. 
The analysis of the upper bound on the transmittance 
has revealed that absorption, radiation, and diffraction 
are three competing mechanisms, responsible for the non- 
monotonic behavior of the transmittance. Furthermore, 
our numerical study has shown that the calculated upper 
bound T n is a good estimate for the actual transmittance 
for all bend radii R > A. In the complementary reversed 
geometry, when the SPPs propagate around a dielectric 
void, it is argued that the fundamental mode is nonra- 
diative, and thus SPP propagation is overall less radia- 
tive. This prediction was verified numerically. Finally, we 
have studied a bend geometry in which a surface plasmon 
resonator is introduced to provide an alternative trans- 
mission channel for enhancing the propagation efficiency. 
Future studies will further analyze resonator design, as 
well as the details of SPP coupling to resonators in the 
presence of loss-compensation. 
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